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— ABSTRfte^-^ 


The inelastic interaction mean free paths X of 5, 10, 

and 15 GeV/c pions have been measured by determining the 

distribution of first-interaction locations in a modular 

tungsten-scintillator ionization spectrometer. In addition 

to commonly used interaction signatures of a few (2-5) particles 

2 

in two or three consecutive modules, we have used a x distri- 
bution to calculate the probability that the first interaction 
occurred at a specific depth in the spectrometer. This latter 
technique seems to be more reliable than use of the simpler cri- 
teria. No significant dependence of X on energy has been observed. 

2 

In tungsten, X for pions is 206 ± 6 g/cm . 
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I . INTRODUCTION 

Studies of hadron-nucleus interactions have long been used 

indirectly to infer information about the fundamental nucleon- 

nucleon interaction. A prime example of these endeavors is the 

consideration of certain subsets of hadron interactions in 

emulsions, such as those with a small number of accompanying 

black tracks, as being associated with the interaction between 

two nucleons. Recently, it has become increasingly evident that 

studies of hadron-nucleus interactions may play a direct role in 

differentiating among proposed models of nucleon-nucleon inter- 
1-3 

actions . 

Apart from its role, direct or indirect, in relation to 
fundamental interactions, the study of hadron-nucleus interactions 
has an important application in reliable energy measurements of 
very high energy particles with ionization spectrometers (calori- 
meters) . In fact, because of the limited past effort in study- 
ing interactions in heavy absorbers, it has been necessary to 
extract information from studies of hadron-nucleon interactions 
in order to construct models of hadron-nucleus interactions, which 
can be applied to the complex physical processes that occur in 
ionization spectrometers . 

Although extensive studies have been made of the probabi- 
lity for interactions of hadrons in hydrogen, relatively little 
information is available on the interaction lengths of hadrons in 
heavy absorbers. This is especially true for tungsten. Accurate 
knowledge of the interaction lengths of nucleons, pions, and 
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heavy nuclei is essential for constructing models of nuclear 
cascades which can be used to extrapolate the calibrated res- 
ponse of ionization spectrometers at accelerator energies to 
the much higher energies to be measured in cosmic rays. 

A tungsten-scintillator ionization spectrometer has been 
exposed to 5, 10, and 15 GeV/c pion and electron beams at the 
Stanford Linear Accelerator Center (SLAC) . Cascade curves for 
pions and electrons, and the energy resolution of the spectrometer 
as a function of incident particle energy, interaction location, 
and spectrometer depth, have been presented in a previous 

7 

publication. We report here results of measurements of the 
inelastic interaction mean free paths of pions at these energies 
in tungsten and discuss techniques for determining the inter- 
action location. 


II. EXPERIMENTAL TECHNIQUE 

The experimental configuration is shown schematically in 

Fig. 1. (See Refs. 7 and 8 for a more detailed description of 

the apparatus and general experimental technique . ) Detectors 

2 

Csl 1 - Cslg were each 2 cm (9 g/cm ) thick CsI(Ti,) crystals. 

Tungsten modules Tg- T^g each contained four layers of tungsten 

2 

(total thickness 79 g/cm ) and three 0.65 cm thick plastic 
scintillators. The ionization was sampled every 3.7 radiation 
lengths (r.l..), but the signal from a module was a measure 
of the ionization over the entire module (11.1 r.l.). Modules 
T l“ T 5' which are callt-d high resolution modules (HRM) , each 
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consisted of 13 g/cm of tungsten followed by a 0.65 cm thick 
plastic scintillator. These modules were used primarily to 
study the early development of cascades, particularly those 
initiated by electrons. For the determination of the pion mean 
free path, signals from , and were added to produce, 

effectively, another thick tungsten module. 

The apparatus was calibrated by determining the response 
of each module to individual non-ir. ter acting 15 GeV/c pions. For 
each event the pulse height of the signal from each module was 
expressed in units of equivalent particles, i.e., the number of 
non-interacting 15 GeV/c pions necessary to produce a signal of 
the same pulse height. The following criteria were used to 
select for analysis only those events in which the beam particle 
did not interact before reaching the tungsten modules: 

a) The response of the upstream scintillator TOF^ was 
required to be between 0.1 and 1.5 equivalent particles. 

b) The response of each Csl module was required to be 
less than 2 equivalent particles. 

Various criteria were then applied to these events to determine 
the distribution of depths corresponding to points of first inter- 
action. The pion mean free path \ was determined from the expo- 
nential tail of this distribution at large depths in the spectro- 
meter. 


HI. INTERACTION CRITERIA 

The first interaction was considered to have occurred at the 
starting point of the nuclear cascade initiated by an incident 
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pion. The parameter measured is the number of the module in 
which the first interaction occurred. Under the assumption that 
the probability for an interaction is independent of energy, the 
probability of an interaction in any module is exponentially 
dependent on the location of the module in the spectrometer. 

A semi-logarithmic plot was made of N(i) vs. i, where N(i) 
is the total number of first interactions occurring in the i*"* 1 
module. A measurement of X was obtained from this distribution 

. q 

using the maximum-liklihood method described by Crannell, et al. 
For each event the probability P for an interaction to occur in 
a particular module is given by the expression 

P = e" t/X (l - e _At/X )/(l - e" L/X ) , (1) 

where At is the module thickness, t is the total thickness up to 
the module in which the interaction occurs, and L is the length of 
the entire spectrometer. The liklihood function £ is defined as 
the product of the individual probabilities P^, determined using 
Eg. (1), for all N events in the data set, i.e., 

X* n p. . (2) 

i*l 1 

Substitution of Eq. (1) into Eq. (2) leads to the relation 

In £ « £ [ + ln(l-e‘ At/X ) - ln(l-e“ L/X )] . 

i«l A 


( 3 ) 
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For each data set in the experiment^ was determined by maxi- 
mizing In £, (which is equivalent to maximizing the liklihood 
function £) as a function of X. Corrections were made for inter- 
actions occurring in the scintillator layers of the modules. 

The difficulty with this method for measuring X lies in 
finding criteria that provide an unambiguous location of the 
first interaction. The interaction of a high energy hadron re- 
sults in the production of secondary particles and a rapid build- 
up .of electromagnetic cascades if neutral pions are produced. 

One usually requires the equivalent of a single particle traversing 
each module upstream of the interaction, plus the equivalent of 
more than one particle in the module in which the interaction 
occurred. The additional requirement that there be more than one 
equivalent particle in at least one module downstream of the inter- 
action reduces the number of erroneous assignments which would 
be made as a result of Landau fluctuations of ionization energy 
loss in a single module. 

The specific signature of an interaction, i.e., the threshold 
number of particles in the interaction module and in subsequent 
modules, selects the physical process to be observed. The use of 
a high threshold, e.g., 10 particles, would not permit detection 
of low multiplicity interactions, unless there were a rapid build- 
up of electromagnetic cascades following the decay of a neutral 
pion into gamma rays. Since higu multiplicity events occur with 
low probability at the energies considered here, it is obvious that 
a very high threshold will result in an excessively large value 
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for the interaction length. On the other hand, the use of a very 

low threshold, e.g., approaching one particle, will result in a 

value for the interaction length which is unquestionably too 

small, because interactions will be simulated by fluctuations in 

the ionization energy losses of single particles. 

Since multiplicities are both material dependent and energy 

dependent, the threshold most suitable for detection of inelastic 

interactions will depend on the material composition of the 

spectrometer modules and on the energy of the interacting hadron. 

The threshold will also depend on the module thickness, since 

this determines the extent of electromagnetic cascade development 

in a single module. We are not aware of any previous measurements 

of the interaction length using a tungsten spectrometer. However, 

several such measurements with energies near the ones we consider 

havti been made with iron spectrometers. Jones, et a_l. ^ have 

reported using a threshold requirement of at least two particles 

2 11 

in two consecutive iron modules 55 g/cm thick. Kurz, et al . 

have used a threshold of at least two particles in three conse- 

2 9 

cutive iron modules 50 g/cm thick. Crannell, et al. have 
tested several criteria and reported best results with the re- 
quirement of three particles in two consecutive iron modules 
2 

50 g/cra thick. 

We have tested several thresholds, including those just 
mentioned. Since the tungsten modules used for our interaction 
length measurements were 11.1 r.l. thick, electromagnetic cas- 
cades reached maximum development in the module in which they 
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originated. Therefore, we obtain better results by using a 
threshold requirement in only two consecutive modules. Specif- 
ically, we have tested the criterion that an interaction occurs 
in the first of two consecutive modules which indicate the 
equivalent of at least two particles. This signature is desig- 
nated as the 2-2 criterion, meaning that there were at least 2 
particles in the interaction module and at least 2 particles in 
the following module. Similarly, we have tested the use of 
3-3, 4-4, and 5-5 criteria for interaction signatures. 

* In addition to these simple interaction signatures, we 
have tried a more elaborate approach in an attempt to improve the 

analysis of the measurements. This approach used the probability 
2 

from a x distribution that an interaction did not occur before 

1.U 

or in the j module. The procedure consisted of the following 
steps: 

1. For each tungsten module j (counting from the upstream 

end of the spectrometer) : 

2 

a) Calculate Xj on the hypothesis that the incident pion 
has not interacted through the j*"* 1 module, using the expression 


j 

E 

i=l 



1 

-) 


2 


( 4 ) 


The symbol NL represents the number of equivalent particles in 
the i th module, and is the standard deviation of the normal 
distribution of signals from single particles traversing that 


module 
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2 

b) From the x distribution with j degrees of freedom, 

2 2 

calculate the probability p^ that x > Xj* i.e., the probability 
that the incident particle has not interacted by the module. 

2. Starting with j = 1, check successive modules until a 

module is found in which the non-interaction probability p^ is 
less than a designated threshold p Q . If 1, disregard this 

module as an interaction module and continue Step 2 for larger 
values of j . If N^> 1, proceed with Step 3. 

3. Check whether the non-interaction probabilities for the 
next three modules are also less than p Q . If so, then module j 
is designated as the interaction module. If not, then disregard 
the module as an interaction module and continue with Step 2 for 
larger values of j . 

The non-interaction probability is required to remain below 
threshold for just three modules following the interaction module, 
rather than for all subsequent modules, in order to avoid biasing 
the interaction location distribution in favor of late inter- 
actions. (A criterion involving all subsequent modules might be 
more easily satisfied for modules farther downstream in the spec- 
trometer, because the non-interaction probability would have to 
remain below threshold for fewer modules.) 

2 

A basic difficulty with the technique is that the Xj are 

2 

not distributed according to the true x “ distribution, because 
the single-particle distributions for each module are not 
Gaussian, as is assumed in the technique. This results in arbi- 
trarily using rather low probability thresholds in order to obtain 
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reasonable results. This problem is not so important for our 
spectrometer, since the added signals from the three scintilla- 
tors in each module tend to make the Landau distribution, 
prominent for a single scintillator, approach a Gaussian distri- 
bution. An empirical test of the technique can be made by 
examining the dependence of the interaction lengths obtained on 
the value chosen for p Q . If the technique has merit, the value 
obtained for the interaction length should be independent of p Q 
over a rather broad range. 

None of the criteria we have used for an interaction signa- 
ture can be expected to indicate for all interactions the correct 
module in which the interaction occurred. Fluctuations in the 
properties of an interaction affect the accuracy of the interaction 
location assignment. Back-scattered particles may result in 
assignments being made in modules upstream of the true location, 
while events with low or zero charged multiplicity may result in 
assignments being made downstream. Arguments have been made to 
show that these competing processes have a negligible effect on 

the relative number of interactions assigned t the central modules 

12 

of a uniform array of identical modules. These arguments are 
valid only if a sufficient number of modules at the two ends of 
the spectrometer are disregarded. The number of modules which 
must be disregarded depends on the material and thickness of the 
modules . We have in all cases disregarded the high resolution 
modules T^- T^ (equivalent to one thick module) at the upstream 
end of the spectrometer as well as the last two modules 
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at the downstream end. This is believed to be sufficient, based 
on the smoothness of the distributions of the number of inter- 
actions occurring in the various modules. 

A typical distribution, obtained for 15 GeV/c pions using 

the x technique with a probability threshold of 0.002, is shown 

in Fig. 2. The number of interactions occurring in each of the 

2 

modules V T 13 was used to determine X. When the x technique 
is used, T^^ is the last possible interaction module, because .-he 
non-interaction probability is required to remain below thr-oi. :ld 
for three modules following the interaction module. (T^g '*‘ 8 
last module of the spectrometer.) 

IV. RESULTS AND DISCUSSION 

The interaction mean free paths X of 5, 10, and 15 GeV/c 
pions in tungsten, obtained using the interaction criteria dis- 
cussed in Sec. Ill, are presented in Table I. For each value the 
quoted error 6X indicates the change in X associated with 
decrease of 0.5 in InjC, from its maximum value. (If li. & is 

9 

normally distributed in X, 6X is one standard deviation. ) If 
the non-interaction probability threshold p Q is high (> .5), 
small fluctuations in the number of equivalent particles in a 
module are identified as interactions. Consequently, a shorter 
mean free path is obtained using such a threshold. On the other 
hand, if p Q is very low, many real interactions are not detected, 
resulting in higher values for X. This increase in X as p Q is 
decreased from 0.001 to 0.0001 can be seen in Table I for pions of 
all three energies. However, the results for each energy are 
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relatively independent of the threshold over a considerable range 
(0.u01<p o ;$0.1) r as required if the technique is to have any 
validity. 

The values of A for 5 and 15 GeV/c pions are consistent 
with each other, but are significantly higher than those for 
10 GeV/c pions. We believe that the lower value,? for 10 GeV/c 
pions may be due to some beam contamination and that it does not 
indicate a real dependence of the mean free path on the incident 
energy. 

For an appreciable number of events the beam particle was 
apparently a muon from a pion decay somewhere upstream of the 
apparatus. If the interaction criteria are not fairly strict, 
these events, because of statistical fluctuations, may occasion- 
ally be designated as interactions. This is equally likely through- 
out the tungsten stack. However, the relati . a number of inter- 
actions which are simulated by muons is greater for the downstream 
modules. Consequently, the tail of the interaction location 
distribution is raised, and a longer mean free path is obtained. 

In order to eliminate events in which an interaction was simulated 
by a muon, a response of three or more equivalent particles was 
required in at least one tungsten module. When this was done, 
the resulting values for A were about 1-2% lower, for probability 
thresholds of 0.05 or less, than the values obtained with no 
discrimination against muons. However, as expected, the differ- 
ences were greater for less stringent interaction criteria 
(higher p Q ) . 
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For the 10 and 15 GeV/c data, the values of A obtained 
using the simpler two-counter criteria exhibit the expected 
threshold dependence, i.e. , the measured values of a increase 
as more particles are required for the interaction signature. 
However, this trend is not observed for the 5 GeV/c data. The 
average multiplicity and, consequently, the number of cascade 
particles and longitudinal extent of a cascade increase with the 
energy of the incident particle. Such interaction signatures as 
the 4-4 and 5-5 criteria are not easily satisfied by cascades ini- 
tiated by 5 GeV/c pions. The requirement of such large rumbers of 
particles, therefore, results in most events being designated as 
non-interacting, and the interaction location distributions so 
obtained are not a good measure of A. 

The fractional error 6 A/A has been used to evaluate the 
various criteria for the determination of A. This quantity, 
obtained in the maximum liklihood calculation, is shown a 
function of p Q if Fig. 3. A well-defined minimum in 6 A/A is seen 
at a probability threshold of about 0.002 for each of the three 
incident pion energies. It can be shown using Eq. (3) that the 

9 

minimum value for 5A/A is given by the expression 

«VXW 2 -° m VV2 ' < 5 > 

It is clear that if the total number of events increases, 6 A/A 
will decrease because of the statistical factor l/»fa. When 
the interaction criterion is fairly loose (high p Q ) , many 
interactions are assigned to the first module, which, in order to 
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reduce boundary effects, is not used in the determination of X. 

As the interaction criterion is made more stringent, fewer inter- 
actions are assigned to the first module and more to the subse- 
quent modules. Consequently, 6X/X in Fig. 3 decreases with 
decreasing p Q , because of the decrease in l/^to. However, when 
the interaction criterion is made too stringent, many real inter- 
actions are missed, so 1//N increases causing 6X/X to increase. 

As the threshold p is decreased from 1.0 to 0.002 (where the 

o 

minima in 6X/X are observed in Fig. 3) , N increases monotonically 
and the l/*fa dependence is no doubt responsible for much of the 
decrease in 6X/X. However, although N continues to increase as 
p Q is decreased further, <5X/X begins to increase. This must indi- 
cate a poorer fit to the hypothesis of exponential decay. Conse- 
quently, we have chosen p Q = 0.002 as the best value for the pro- 
bability threshold. 

For the 5 and 10 GeV/c data, the 6X/X for the two-module cri- 
teria are considerably larger than the minimum values obtained 
2 

using the x technique. For the 15 GeV/c data, the 6X/X are com- 
parable or slightly larger for the two-module criteria. Conse- 

2 

quently, for the energies considered in this experiment, the x 
technique seems to be the best method for determining X. 

The mean free path for pions in tungsten and other heavy 
absorbers can, of course, be determined by methods which do not 
employ an ionization spectrometer. The absorption cross section 
°abs can be measure<a b y standard techniques and X calculated 
according to the relation 
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X 


N 0 ^ABS 


( 6 ) 


where A is the atomic mass number of the absorber and N q is 
Avogadro’s number. For the purpose of comparing data for dif- 
ferent absorbers, it is more informative to consider O^g rather 
than X as a function of A. According to a simple opaque sphere or 
"billard ball" model, CT^gg should be proportional to A ' , if 

the target nucleus is perfectly opaque to the incident particle. 

However, absorber nuclei are to some extent transparent, more so 

3 /4 

for lower A, resulting ir.. a steeper A dependence (more like A ' ) . 

2 

The mean free paths obtained using the x technique with 

p Q = 0.002 are nearly the same for 5 and 15 GeV/c pions. If these 

2 

data are combined, our result is a value of 206 + 6 g/cm for 

the pion mean free path in tungsten, or equivalently, a cross sec- 

13 

tion a = 1480 + 45 mb. Allaby, et al. have determined the absorp- 
tion cross sections for a series of absorbers, using the method of 

removal of particles from an accelerator beam. A value of 

2 

X = 214 + 3 g/cm (a = 1425 + 20 mb) is obtained for tungsten by 
interpolating between their values (20 GeV/c it ) for tin and lead. 

A value of X = 207 + 6 g/cm^ (o = 1475 + 40 rob) is obtained for 

tungsten by interpolation between the values obtained by Denisov, 

14 9 

et al . for 13.3 GeV/c pions in tin ind lead. Crannell, et al . 

2 

quote a value of X = 164 + 9 g/cm (O' = 565 + 31 mb) for 17.8 
GeV pions in iron. Their experiment involved an ionization 
spectrometer somewhat similar to ours. 

Measured pion cross sections obtained in these and other 
experiments^ are shown as a function of A in Fig. 4. Our result 
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for tungsten is consistent with the values obtained for different 
absorbers in the other experiments. The straight line shown in 
Fig. 4 was determined by a non-weighted least squares fit to all 
the data. It can be seen that the data fit the hypothesis 
a — aA~ # where a = 31 miiiibarns and b = 0.74. 

Y. SUMMARY AND CONCLUSIONS 

The interaction mean free path X of pions in tungsten has 
been determined from the exponential distribution of starting 
points of cascades in a tungsten scintillator ionization spectro- 
meter. The value obtained for X depends on the specific signa- 
ture used to locate the depth of the first interaction from which 
the cascade ensues. If the interaction signature is based upon a 
low number of particles, fluctuations in the energy loss of single 
particles may falsely indicate an interaction. If the signature 
is based upon a large number of particles, cascades initiated by 
low energy primaries may not contain enough particles to satisfy 

the requisite signature, and some interactions will not be detected. 

2 

The technique utilizing the probability from a x distribution to 
locate the depths in the spectrometer of the first interactions 
seems to provide a less biased measurement of X, since a rather 
wide range of probability thresholds gives essentially the same 
value for X. By selecting the probability threshold which mini- 
mizes the fractional error 6X/X , we have obtained a value of 

2 

X - 206 ± 6 g/cm as the best estimate for the inelastic inter- 
action mean free path of pions in tungsten. This value is in 

agreement with the absorber dependence of the interaction cross 

3/4 

section given by a ■ 31 A ' 


miiiibarns. 
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TABLE I 

Mean free paths X of pions in tungsten. The probability thresholds 

2 

p Q which were used in the x technique and the two-module criteria, 
as defined in the text, are given in the first column. 


. 2 
P 0 for x 

Technique 


0.5 

0.2 

0.1 

0.05 

0.02 

0.01 

0.005 

0.002 

0.001 

0.0005 

0.0002 

0.0001 

Two-module 

Criteria 


2-2 

3- 3 

4- 4 


X (g/cm 2 ) 


5 GeV/c 

10 GeV/c 

15 GeV/c 

187 ± 

16 

176 

± 

10 

201 

± 8 

192 ± 

11 

187 

± 

11 

204 

± 7 

198 ± 

12 

188 

± 

11 

205 

± 7 

203 ± 

12 

186 

± 

10 

206 

± 7 

202 ± 

12 

186 

± 

10 

207 

± 7 

204 ± 

12 

188 

± 

10 

206 

± 7 

205 ± 

12 

187 

± 

10 

206 

± 7 

207 ± 

12 

190 

± 

10 

206 

± 7 

209 t 

12 

192 

± 

10 

207 

± 7 

211 ± 

12 

195 

± 

11 

209 

± 7 

215 t 

13 

199 

± 

11 

214 

± 7 

217 ± 

13 

203 

± 

11 

216 

± 7 


198 

± 

12 

205 

± 

12 

219 

± 

8 

201 

t 

13 

217 

i 

12 

228 

± 

8 

198 

1 

14 

225 

± 

13 

226 

± 

8 

188 

± 

16 

224 

t 

14 

232 

± 

8 


5-5 
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FIGURE CAPTIONS • 

1. Scale drawing of apparatus. SC^, SC 2 and SC^ are wire spark 
chambers; TOF^, T0F 2 ' B l' P 2' an< * P 5 arF pl ast i c scintillators. 

2. Frequency of first interactions at various depths t in the 

_ 2 

spectrometer for 15 GeV/c pions. The x technique was used 
with a probability threshold p o c 0.002. The errors are 
statistical (cr//N) . 

3. Fractional uncertainty 6X/X in the mean free path of pions in 

2 

tungsten. The x technique was used with different probability 
thresholds p Q . Results at 5 GeV/c; a. Results at 10 GeV/c: 

O . Results at 15 GeV/c: A. 

4. Absorption cross section <J vs^. atomic mass number A of the 

absorber. Results from this experiment: • . Results from 

Ref. 9: ■ . Results from Ref. 13: A . Results froir Ref. 14: 

V . Results from Ref 15: 0 . The straight line was deter- 

mined by a non-weighted least squares fit to all the dtta. 
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